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Dewetting of thin polymer films near the glass transition 
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Dewetting of ultra-thin polymer films near the glass transition exhibits unexpected front mor- 
phologies [G. Reiter, Phys. Rev. Lett., 87, 186101 (2001)]. We present here the first theoretical 
attempt to understand these features, focusing on the shear-thinning behaviour of these films. We 
analyse the profile of the dewetting film, and characterize the time evolution of the dry region ra- 
dius, Rd{t), and of the rim height, hm{t). After a transient time depending on the initial thickness, 
hm{t) grows like \^ while Rd{t) increases like exp(\/?). Different regimes of growth are expected, 
depending on the initial film thickness and experimental time range. 

PACS numbers: 68.60.-p, 68.15.+e, 68.55.-a, 83.10.-y 



Thin liquids films are ubiquitous entities in science and 
technology. In engineering, for instance, they serve to 
protect surfaces, and applications arise in paints, ad- 
hesives and membranes [l|. Thin hquids films display 
a variety of interesting dynamics phenomena and have 
therefore been the focus of many experimental and the- 
oretical studies ||]. Nearly half a century ago, Taylor ||] 
and Culick |^ analyzed the growth of a circular hole in a 
thin liquid sheet Q . By balancing surface tension forces 
against inertia, they found that the rim of liquid at the 
edge of the films retracts at a constant velocity, a predic- 
tion first checked experimentally by Mc Entee and Mysels 
1^. The precise shape of the rim was later analysed by 
Keller et al. [Q, ^ who showed that it is a cylindrical 
cap expanding in time like ^/i. Recent experiments by 
Debregeas and collaborators on thin suspended films 
of very viscous liquids revealed unexpected features: (a) 
First, the retraction velocity grows exponentially with 
time (with a characteristic time scale = hiri/\S\, where 
hi, Tj and S are respectively the initial film thickness, the 
viscosity and the spreading coefficient (b) second, 

the liquid is not collected into a rim and the film remains 
flat through the retraction. According to these authors, 
the uniform thickening of the retracting film was a con- 
sequence of its viscoelasticity, which permits an elastic 
propagation into the film of the surface tension forces act- 
ing on the edge. Solving numerically the Navier-Stokes 
equations for long wavelength modulations of the film, 
Brenner and Gueyffier showed, however, that the ab- 
sence of rim can also result from a purely viscous effect. 
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FIG. 1: Film geometry: h(r,t) is the profile of the film, hmit) 
is the height of the rim, and Rdit) is the radius of the dry 
zone. The initial {t = 0) step-like profile is represented by the 
dashed line. v(r,t) is the radial, axisymmetric fiow field [p^ . 



a much slower decay on the rear side, builds up progres- 
sively In this Letter we present the first theoreti- 
cal attempt to understand these new features. Figure |l| 
shows the film geometry. 

In order to characterise the rheologic properties of the 
film, we introduce the stress tensor cr^ = ~p5ij + (t™^ , 
where p is the pressure and cr™y represents the effects of 
internal friction. We relate fi™ ,- to the strain rate tensor 
7y by a constitutive law of the form: 



'(7) = f^o«'(7T) 



(1) 



Very recently |1^, Reiter studied the dewetting of ul- 
trathin (i.e. thinner than the coil size), almost glassy 
polystyrene (PS) films deposited onto silicon wafers 
coated with a polydimethylsiloxane (PDMS) monolayer. 
He found that a highly asymmetric rim, with an ex- 
tremely steep side towards the interior of the hole and 



where ctq and t are material constants, and $ is a 
generic function. For a purely viscous liquid, the function 
<i> is linear, whereas for polymers, just above Tg, it can 
be shown within the framework of the free- volume model 
|l5i that cr'" is expected to vary logarithmically with 7 
ailll: 



$(7r) = ln(l + 7t) 



(2) 
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At low strain-rates (7r < 1), this law displays a 
viscous- type behaviour (cr™ « 7707, with a zero-shear vis- 
cosity 770 = o-qt), while for large values of 7T, cr™ reaches 
an almost constant value (shear-thinning behaviour). 
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FIG. 2: Film profile for K = O.l/i* and t = 10" V. For clarity 
reasons, a sectional view is displayed on a quarter of circle. 



Assuming the fluid to be incompressible, mass conser- 
vation leads to: 



1 ^dh{r,t) 



rf3{r,t) dh{r,t) a{r,t) - f3{r,t) 



h{r, t) ' dt 



dr 



(3) 



Equation ^ involves two unknown functions a{r^ t) 
and I3{r,t)^ that are positive, dimensionless forms of the 
strain-rate components 7rr and : 



dr 



(4) 



We thus need two additional equations to determine 
h{r,t). First, note that the following partial differential 
equation can be directly derived from Eq.[|: 



dr 



P + a 



(5) 



Neglecting the inertial term, conservation of momen- 
tum (projected on the radial direction leads to: 



dar 



dr 



= 



(6) 



It can be shown that Eqs. |^ and ^, along with the free- 
surface boundary condition (i.e. a^z = at the contact 
with ambient atmosphere allow one to express the 

strain rate a as a function of /3 only : a = F{P). Substi- 
tution of the constitutive law Q in Eq.^ then leads to the 
following differential equation for F: 



dp 



F{p) + P 



(7) 
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FIG. 3: Evolution of the film profile h{r,t) for hi — O.lh* 
and t ranging from t = to t = 5.10~^r with a time step of 
10"V. 



which has to be solved along with the condition F{0) = 
(far away from the perturbed central region, the strain 
rates must decrease to zero). 

In order to solve Eq.^, we should supply it with initial 
and boundary conditions. Our initial profile is assumed 
to be uniform, with a thickness hi, except in a bored 
region ranging from r = to r = Rq (as shown in fig.^. 
Rq is our characteristic radial length used thereafter to 
make r dimensionless: 



h{r,t^ 0) 



hi 




if r > i?o = 1 
otherwise 



(8) 



We do not consider here the origin of the initial dewet- 
ting process: experimentally, it is found that a thick PS 
film on a silicon substrate is metastable and dewets via 
nucleation and growth of dry patches [ p^ , while thinner 
films {hi <100 nm) are unstable and dewet by spinodal 
decomposition j2^. In his latest experiments on ultra- 
thin films [|l2j, Reiter characterized the early stage of 
the dewetting process by the formation and coalescence 
of little holes, with the displaced material uniformly dis- 
tributed between the holes, without visible rims. Our 
initial time t = might correspond to the end of this 
preliminary process. 

Equation ^ applies outside the dewetted region 
(Rdit) < r < 00). At the edge of the rim (r = Rd{t)), 
the rim height, hm(t), can be determined by taking into 
account capillary forces. The applied force on the rim, 
pushing the film away the dry area, must be balanced 
by the internal radial stress: |5'| = |crrr|^m(i). Assum- 
ing the lateral extension of the film to be large enough, 
the film thickness must reach hi far from the dry region: 
limr^oo h(r, t) — hi (Vi). The complete resolution of our 



set of equations 
characteristics p5|. 




can be achieved using a method of 
Thereafter, all thicknesses will be 
made dimensionless by normalizing with a characteristic 
length h* = \S\/aQ (for Reiter's experimental conditions, 
we estimate h* to be of the order of 500 A). 

Note first that for a purely viscous liquid {^{jr) = jr), 
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FIG. 4: Logarithmic plot of the dry zone radius Rd versus 
\/t/T, for hi — 0.1, 0.5, 1, 5 and 10 h* . Note the hnear 
behaviour of log [Rd{t) / Rq\ with \/t/T after a transient time. 



our model leads to a constant and uniform thickness for 
the film, with an exponential growth of the dry radius: 



Rd(t) = Roe' 



Roe' 



(9) 



This is in complete agreement with the experimental 
results of Debregeas et al. for the dewetting of suspended 
polymer films |^ and supported films plj l. 

If we take into account the shear-thinning behavior of 
the polymer film (o'™(7) = (Tq ln(l + -yr)) the dewetting 
front is quite different. An example of profile obtained 
for hi = O.l/i* is shown in Fig.^. This profile is charac- 
terized by a highly asymmetric shape for the rim. There 
is a striking similarity between such a profile and those 
observed by Reiter for ultra-thin films |l^. On fig.^ we 
present the time evolution of this profile (for t ranging 
from t = Q to t = 5.10~^t, with a time step of 10~^t). 
The velocity field v{r,t) is found to be almost propor- 
tional to 1/r at any time for thick profiles {hi ^ /i*), 
similar to the long range radial plug flow observed by 
Debregeas et al. We now focus on two natural pa- 
rameters characterizing the film profile: the dry radius 
Rd{t) and rim height hm{t). 

In order to discuss the time dependence of the rim 
height, we need to compare hi with the characteristic 
thickness h* (in Fig. |^ we present our results for differ- 
ent initial thicknesses). For very thin profiles {hi h*), 
the increase of rim height is initially very fast, whereas for 
thicker films {hi ^ h*), hm{t) is nearly constant during 
a long period of time, before growing faster. For the sake 
of conciseness, we will mainly focus on the case of thick 
films, for which simple scaling laws can be analytically 
obtained. As shown in Fig.||, for hi — 100ft,* two differ- 
ent time regimes can be distinguished: at the beginning 
of hole formation, h^ is nearly constant and equals its 
initial value hi. We have checked numerically (cf. Fig.^ 
b) that the anticipated analytical behaviour 



h* 



is obeyed after a crossover time to ~ '''{j^)'^ 
In Fig.^, we present the time evolution of Rd for differ- 
ent initial thicknesses hi . It is physically understandable 
that Rd, at a given time t, is larger for thinner films : for 
a given applied force \S\ per unit length, the thicker the 
film, the more the material to be displaced and the lower 
the dewetting velocity. Fig.^ shows that the following 
anticipated analytical behaviour is obeyed for to: 



Rd{t) ~ Ro 

t»to 



(11) 



Note that the coefficient R^^ strongly depends on hi. 
During the early stage {t <^ t^), the strain rates at the 
rim are constant and small, as a consequence of the con- 
stant large thickness. In this case, the rheological law 
of the film is viscous-type and the dry radius increases 
exponentially, as for the case previously discussed (see 
Eq.^. The range of thicknesses experimentally studied 
by Debregeas et al. (5 to 250 /im) covers the domain 
hi ^ h* . As the corresponding crossover time is very 
large (tg 3> t), these experiments only covered mainly 
the first regime t <ti tQ (this criterion hi ^ h* is in fact a 
basic hypothesis of their "soft balloon" model |^). The 
crossover to the exp(-\/t) regime is a consequence of the 
non-linearity of rheological law ^j. 

The case hi <C h* will not be fully discussed here (cf. 
JlSf). Briefly, we can note on Fig.^that a sharp rim grows 
very quickly (with an initial speed hm ~ hiSxph* /2hi) 
at short times. In this regime, hm{t) is well described by 
the integral equation: 



-da; 



t<r r 



(12) 



(10) 



Later on, h^ and Rd reach the same laws of growth as 
in the case of thick films (Eqs.^ and |l^). 

In conclusion, our model accounts well for the expo- 
nential growth and the absence of rim characteristics of 
the dewetting regime observed by Debregeas et al. for 
viscous polymer films. Taking into account the shear- 
thinning behaviour of the polymers near Tg enables us to 
see the modifications induced by this particular rheology 
on the film morphology |^ . It appears from our results 
that the early stages of dewetting for a shear-thinning 
polymer film are mainly determined by the ratio of its 
initial thickness hi to a characteristic scale h* (related to 
the driving force of the process, S, and the rheological re- 
sponse of the material, characterized by ctq)- In all cases, 
after a transient time depending on hi, the rim height 
grows proportionally to \/t/T while the dry radius ex- 
pands like exp (4-\/t/r). The film profile exhibits a sharp, 
asymmetric rim similar to the one observed by Reiter in 
his latest experiments with ultrathin PS films. Work is 
now in progress to incorporate in our model the Laplace 
pressure (se e ]l8[) which might lead to oscillatory dewet- 
ting fronts |11|, |13[, and to evaluate the possible effects 
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FIG. 5: Time evolution of the rim height hm'- (a) Plot of hm versus t/r, for hi = 0.01, 0.1, 0.5 and 1 h* . Note that the curves 
corresponding to hi =0.1 and O.Ol/i* cannot be distinguished within the graph scale, (b) Plot versus \pijr to show the linear 

behaviour of hm with \Jt/T after the transient time to- For a thick film (hi ^ h*), two regimes can be distinguished. For 
i < to ~ ^ij^Y , hm is constant and Rd{t) increases exponentially with time. For t > to, hm growths linearly with ^/Tp^ and 
Rd{t) expands like exp \ftjT. 



of chain entanglements. In the future, we aim to study the material may lead to new dewetting morphologies, 
the dewetting of thin polymer films below Tg, where the We wish to acknowledge fruitful discussions with A. 
existence of a yield stress in the rheological response of Arabian, G. Debregeas, and G. Reiter. 
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